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One-dimensional transistors are well-understood today; computer 
techniques for detailed large-signal and small-signal analyses are available. 
Real transistors are three-dimensional, however, and lateral effects are 
only understood qualitatively. Accurate modeling of lateral effects cannot 
be accomplished without quantitative analyses of three-dimensional trans- 
istors. Unfortunately, even the simplest analysis of lateral effects leads to a 
partial differential equation. In this paper, a fast and accurate numerical 
technique is used to solve the partial differential equation. This makes 
feasible a three-dimensional small-signal analysis of transistors operating 
in the low-injection regime. 

Calculated h-parameters for a high-frequency, double-diffused, silicon 
transistor are in good agreement with experimental values. 

I. INTRODUCTION 

One-dimensional transistors are well-understood today; computer 
techniques for detailed large-signal and small-signal analyses are 
available. 1 A new charge-control model, 2 which is quite promising for 
use in circuit analysis programs, has been developed with the aid of 
insight obtained from the large-scale computer analyses. 

However, real transistors are three-dimensional and lateral effects 
are only understood qualitatively. Accurate modeling of lateral effects 
cannot be accomplished without quantitative analyses of three-dimen- 
sional transistors. Unfortunately, even the simplest analysis of lateral 
effects leads to a partial differential equation [namely (7)]. Except for 
certain very simple transistor geometrical configurations, an analytic 
solution to the partial differential equation is not possible, and up to 
now, numerical solutions have been at best difficult and expensive in 
computer time. 

In this paper, a fast and accurate numerical technique 3 ' 4 is used to 
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solve the partial differential equation. This makes feasible a three- 
dimensional small-signal analysis of transistors operating in the low- 
injection regime. 

It is hoped that with the aid of this technique, quantitative under- 
standing of lateral effects in transistors can be achieved, and an accurate 
but simple model can be developed. 

In Section II, the partial differential equation for the potential in the 
base region of the transistor is derived. Section III briefly reviews the 
solution technique. In Section IV, the calculation of /t-parameters from 
the partial differential equation solution is described. Section V presents 
the results from a sample calculation, and compares them with experi- 
mental values. 

II. TRANSISTOR MODEL 

In this section, a discrete bipolar transistor with collector contact 
made to the substrate is analyzed. (A planar bipolar transistor could 
equally well have been chosen.) The small-signal behavior of the 
transistor in one dimension is modeled and a partial differential equation 
in the other two dimensions is derived. 

Figure 1 is a cross section perpendicular to the surface of a typical 
discrete transistor with diffused base and emitter. A partial differential 
equation for <p, the potential in the base region, will now be derived. 
Let z be the coordinate perpendicular to the surface of the transistor. 
Let <p(x, y) be the potential in the base, p(x, y) be the charge density, 
and J(x, y) be the current density parallel to the transistor surface; 
<p, p, a, and J are average values, and are taken to be independent of z. 

Ohm's law is 

V<p=-~J, (1) 

<r 

where a is the conductivity, and the continuity equation for the charge is 
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Fig. 1 — Cross section of transistor. 
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v-J + ff-o. (2) 

The dp/dt term is composed of capacitatively and resistively injected 
charge into the base, and will be considered shortly. Taking the diverg- 
ence of (1) and substituting (2), one obtains 

2 d*tp d 2 tp 1 dp . 

Bx By a Bt 

In order to calculate Bp/Bl, the net rate of the charge injection into 
the base region, a specific model is necessary. This paper will consider 
only a simple model, since the purpose is to illustrate the technique 
rather than to model exactly some particular transistor. Two separate 
sets of simplifying assumptions will be made, one set in calculating 
I/a Bp/Bt to set up the partial differential equation, and the second set 
in calculating small-signal parameters from the solution of the partial 
differential equation. The second set will be discussed in Section IV. 

It is assumed that the emitter doping is sufficiently high that the 
potential on the emitter side of the emitter-base junction is constant. 
This constant is taken equal to zero. It is similarly assumed that the 
potential on the collector side of the collector-base junction is constant. 
The base region is divided into two parts — the active base region, directly 
under the emitter diffusion, and the passive base region, the remainder 
of the base. The conductivity of each of the two regions is taken to be 
constant, but the two conductivities are not equal. Finally, low injection 
is assumed so that current and charge injected by the emitter-base 
junction is uniform in the active base region. 

With these assumptions, the rate of charge injection into the base may 
be calculated. Equations (4), (5), and (6) may be interpreted with the 
aid of the equivalent circuit of Fig. 2, which is essentially a hybrid-ir 
model of a transistor neglecting base resistance corrections. The charge 
injected into the base at point (/, y) may be obtained by calculating 
the base current for the circuit in Fig. 2; such an equivalent circuit 
obtains at each point of the active base. Only the circuit elements which 
contribute to charge flow in the z-direction are included; the lateral 
charge flow is calculated by solving the partial differential equation (3). 
Thus, it is not necessary to include in the circuit of Fig. 2 the usual 
hybrid-7r base resistances. In the equations below, ip(x, y) = v BE (x, y) 
is the only quantity that depends on x and y, for (j*, y) in the active 
base region. The rate of charge injection into the base has three compo- 
nents. 
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Fig. 2 — Equivalent circuit used to calculate charge injection into the base and 
current injection into the collector, at a point (x, y) of the active base region. 



(t) The dynamic resistance and capacitance of the emitter-base 
junction contributes to dp/dt only in the active base region. 



8). -A (*£+*')' 



(4) 



where W A is the width (z-direction) of the active base region; 
C E is the dynamic capacitance per unit area of the emitter-base 
junction, and includes the ordinary junction capacitance; and 
G E is the conductance per unit area. 
(it) The collector junction capacitance contributes 



§ at)=wit {(p - Vo * ) > 



(5) 



where C c is the capacitance per unit area of the collector-base 
junction, v CE is the potential of the collector side of the collector- 
base junction, and W in the active base region is W A and in the 
passive base region is W P , the width of the passive base. 
(Hi) Base-width modulation in the active base region contributes 



(ti. - £(*+*£)* -** 



(6) 



where the conductance G' and capacitance C are each per unit 
area. Generally C and G' are not important. In addition, base- 
width modulation is the source of G , a conductance per unit 
area. 

In the passive base region, the only charge injected is through the 
capacitance C c • 

Collecting the three terms, assuming frequency dependence e' ut , 
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defining sheet resistance R A = l/W A <r A {mdR P = \/W P a P , the following 
partial differential equations for the potential are obtained. 

Active base: 

V> = R A [G E + G' + u»(C c + C E + C')\<p 

- R A [G' + ua(C + Cc)]vce ■ (7) 
Passive base: 

VV - io>RrC c (<p - v CE ). (8) 

For typical parameters, 

G E » Go, 

C, » C c + C, (9) 

C c »C"; 

then in the active base the equation simplifies to 

VV - Ra[(G e + ««C*) V - ((?' + iuC )voml (10) 

This equation was used for numerical calculations. 
Equations (8) and (10) are of the form 

VV = 7 2 (v - Vce), (11) 

where y 2 is frequency-dependent and takes on different values y A and yl 
in the active and passive base regions, but otherwise is independent 
of x and y. 

Any other modeling of the z-direction of the transistor is equally 
usable, provided that it leads to equations of the type 

VV = a<p + b, (12) 

where a and b are independent of x and y in each of several regions. 
Figure 3 shows a typical double-base-stripe discrete bipolar transistor 
from the top. The active base region is under the emitter diffusion, 
marked E; this region is represented by eq. (10). Under the two regions 
marked B, ip is constant and <p = v BE . The remainder of the region 
shown is the passive base region; in this region, eq. (8) holds. Around 
the passive base region is the boundary of the base diffusion. A boundary 
condition is that no current flows through this boundary; the current 
flowing through the boundary is proportional to the normally-directed 
derivative of the potential, so on this boundary 

n-W = |; = 0, (13) 
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Fig. 3 — Top view of transistor used in calculations. 

where n is the outer-directed unit normal to the boundary. At the edge 
of the emitter diffusion, the dividing line between the active and passive 
base regions, the boundary conditions are that the potential be contin- 
uous and that the total current flow be continuous. These conditions are 



fA = <Pr; 



n (i: v ^) = n fe v ^)' 



(14) 
(15) 



where n is the unit normal outward from the emitter diffusion, <p A refers 
to the potential just inside the active base, and tp P refers to the potential 
just inside the passive base. Equation (14) neglects the edge injection 
of current from the emitter-base junction; such injection could easily 
be included. It leads to a boundary condition, similar to eq. (15), but 
involving both p and dip/dn. 



III. SOLUTION TECHNIQUE 

This section describes an economical method for solving the coupled 
partial differential equations (8) and (10) with boundary conditions (13) 
through (15). The usual way to solve such equations is by finite dif- 
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ferences. A square or rectangular grid is superimposed on the x, y plane 
and derivatives at the grid vertices are approximated by finite differences. 
A large number of linear algebraic equations are generated and solved 
iteratively. While very powerful, finite difference methods have several 
disadvantages for the problem at hand. 

It may not be easy to get the transistor boundaries (emitter diffusion, 
base contacts, and base diffusion) to fall on grid lines, so that more 
complicated programming is necessary. The areas in which the equations 
must be solved are frequently odd-shaped (see Fig. 4). In such cases, 
the usual iteration schemes are not guaranteed to converge and it may 
be necessary to develop an efficient scheme. 

The method 3,4 developed for solving equations like (11) does not 
employ a grid and does not require iteration to solve the set of linear 
algebraic equations. The method will be briefly described for a single 
equation, 

VV = yVi (16) 

holding in some area A with boundary T. The area A may be any shape; 
it may even be multiply connected. For (16), Ref. 4 derives the exact 
expression, 



r<p(s') 



= F l[ K ° 



( ' r) ^ + (n " r/r)Kl(7rMs) ] 



(Is. 



(17) 



Here s' is any point on the boundary (except at a corner point), P 
denotes the Cauchy principle value, K n and K\ are modified Bessel 
functions of the first kind, 5 and r is the vector from s' to s. Points inside 
the region A are not involved; however, if s' is inside A, (17) can be 



Fig. 4 — "Exploded" view of one-fourth of transistor. Left — passive base region; 
right — active base region (under emitter). Numbers indicate boundary conditions 
applying (see text). 
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used if v is replaced by 2t. Thus, if <p and dip/dn are known on the 
boundary, cp may be found inside the boundary. 

To find <p(s) and d<p(s)/dn on the boundary, (17) is regarded as an 
integral equation. This and the boundary conditions for the original 
partial differential equation suffice to determine the solution. In prin- 
ciple, an accurate numerical solution is quite simple, although the 
computer programming is rather complex. Various methods may be 
used to solve (17) plus the boundary conditions; the method used in 
Ref. 4 will be described here. The boundary T is divided into a number 
of straight-line segments; on each segment <p(s) and d<p(s)/dn are approx- 
imated by polynomials of any desired order with unknown coefficients. 
Using geometrical information about r, and the expansion of K a and Kj 
in terms of logarithms and polynomials, 5 the integrals in (17) may be 
done symbolically. 4 Equation (17) cannot in general be satisfied exactly 
at each point s' along the boundary; N points are chosen at which to 
satisfy (17) exactly, where N is the total number of unknown coefficients 
in all of the polynomials approximating <p and d<p/dn, on every part of r. 
Then (17) is replaced by N linear algebraic equations for the poly- 
nomial coefficients. For typical transistors, N is no more than 50 or so, 
and therefore the equations may be solved directly without iteration 
by Gaussian elimination; N = 50 gives an accuracy of about one percent 
for the numerical example considered later. 

When the polynomial coefficients have been found, ip and d<p/dn are 
known on the boundary and <p may be found inside A if it is desired. 
For the transistor analysis, however, <p is not needed inside A; it is only 
needed on the boundary. Although the program is quite complex, the 
computer time is quite short. For Laplace's equation, the boundary 
integral equation method is considerably faster than finite difference 
methods. 4 

The two coupled partial differential equations, (8) and (10), may be 
solved in a similar manner. Fig. 4 shows the areas of Fig. 3 in which 
each of the two equations can be used. Because of the fourfold symmetry 
of Fig. 3, only one-fourth of the transistor need be considered, and only 
one-fourth is shown. On the left is the passive base region, in which (10) 
holds; on the right is the active base region, in which (8) holds. The 
numbers on the line segments indicate the boundary conditions on each 
segment. On segments numbered 1, the base contacts, (<p — v C s) is 
constant; on the exterior of the base diffusion, segments numbered 2, 
and symmetry lines, dip/dn = 0; and on segments numbered 3, the 
boundary of the passive and active base regions, (14) and (15) hold. 
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IV. CALCULATION OF /l-PARAMETERS 

With a fast and accurate method of solving the coupled partial 
differential equations (8) and (10), finding the small-signal parameters 
of the transistor is relatively easy. This section describes the method 
for determining common-emitter /(-parameters. 

Common-emitter /t-parameters are defined by 

v B e = hui B + h 12 v CE = h ie i B + h re v CE , (18a) 

ic = h 2l i B + h 22 v C E = h fe i B + h oe v CB . (18b) 

In order to find all four /(-parameters at a given frequency, two 
potential solutions are necessary. The first solution gives h n and h 2l 
and assumes v CE to be zero. With v CE = 0, v BE = h u i B and i c = h 2X i B ■ 
The base contacts are taken to be at unit potential, v BE = 1. The bound- 
ary integral equation method is used to solve (8) and (10) with boundary 
conditions (13), (14), (15), plus <p = 1 on the base contact (lines num- 
bered 1 in Fig. 4). In Fig. 4, the output is d<p/dn on lines numbered 1, 
where ip = 1; <p on lines numbered 2, where dip/dn = 0; and both <p and 
d<p/dn on lines numbered 3. The current flowing into the base contact 
is just 



In 



-11^"^ 09 



where the integration is over lines numbered 1 in Fig. 4. Since d<p/dn 
is approximated by a polynomial on each of the sides, the integral 
in (19) is trivial. This gives h u = \/i B . 

To find h 2 \ , the collector current must be found, which involves the 
modeling of the transistor in the z-direction. Again, the simple model 
used may be interpreted with the aid of Fig. 2. Transit-time effects in 
base and collector will be neglected. The z-directed collector current 
flowing through a small area Ax Ay has three components: 

(i) Diffusion of injected charge across the active base region. 
G m is a conductance per unit area. 

(A« c )i = G m <p Ax Ay. (20) 

(ii) Injected current from collector-base junction capacitance. 

(A* e )a = — C c ^ (<p — v CE ) Ax Ay. (21) 

(in) Injected current from base-width modulation in the active base 
region. 
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(A» c ). = GoVcb Ax Ay - ((?' + C |)(*> - v CE ) Ax Ay. (22) 

From the passive base region, current flows to the collector only 

through the capacitance C c ■ 
The total collector current is 

i c = [f w« - <n - MCc + c)k 

J J Active Bnse 

+ [(Go + G') + MGc 4- C')]v C e} dx dy 

— fl iwc(<P — v CE ) dx dy. (23) 

•'•'PiisBive Base 

In the calculations, G' was neglected, G + G' = G t were combined, 
and C was neglected. 

Thus it is necessary to do an area integral of the potential. Since the 
factors multiplying <p are independent of x and y, a two-dimensional 
numerical integration may be avoided. Consider one of Green's bound- 
ary-value formulas, 

Jl «?W - *V'G) dx dy = f T (a | - , f n ) *, (24) 

where the line integral is over T, the boundary of the area A. This is 
true for any well-behaved functions <p(x, y) and G(x, y). Suppose that 
<f> obeys (16). In the area integral, substitute y 2 <p for VV to obtain 

-11 *7>0 - y'G) **-£(<»£-* f) *■ (25) 

Now G may be chosen so that 

V 2 G - 7 2 G = 1, ( 26 ) 

or 

G = -lM 
to give 

/J>* *--?££>■ (27) 

Thus the area integral in (23) may be replaced by a line integral; 
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since d<p/dn is a known polynomial, the integral is again trivial. In this 
way i c may be found, and h 2l = icI^B ■ 

In order to find h l2 and h 22 , Ib is assumed to be zero. With i B = 0, 
v B e = h 12 v C E and i c = h 22 v CE . The potential at the base contacts, v BE , 
is unknown; since i B = 0, d<p/dn = at the base. For the second solution, 
the boundary integral equation method is used to solve (8) and (10) for 
(<p — v CE ) with boundary conditions (13), (14), (15), plus <p = v BE = 
unknown constant and d<p/dn = on the base contact (lines numbered 1 
in Fig. 4). v CE is assumed to be 1. The output is <p on lines numbered 1 
and 2, where d<p/dn = 0, and both <p and d<p/dn on lines numbered 3. 
Then v BE = h l2 '\s found directly from <p on the base contact, and i c may 
be found as before to give h 22 = i c . 

V. RESULTS FROM SAMPLE CALCULATION 

Experimental /i-parameters were available for a double-diffused, 
double-base-stripe silicon transistor, similar to that shown in Fig. 3. 
The data were taken with v CE = 1 V, and I E = 10 mA (approximately 
230 A/cm 2 ). In the theoretical calculations, there are seven parameters: 
R P , C c , Ra , Gm , G E , C B , and G t . Two of these, R P and G m , were not 
adjusted. R P , the sheet resistance of the passive base, was given its 
nominal fabrication value. G m , which is essentially a transconductance 
per unit area, was set equal to I E /(A E kT/q). A E is the emitter area. 
The other five parameters were chosen so that the experimental and 
calculated /i-parameter matched at low frequencies (below 1 MHz). 
No additional adjustments were made to fit the high-frequency data. 
The parameter values used in the calculation are given in Table I. For 
the calculations, the nominal fabrication geometry was used. 

Complex-plane plots of the four common-emitter ^-parameters are 
given in Figs. 5 through 8. The solid lines are drawn from the experi- 
mental data, only a few points of which are shown (solid circles). Fre- 



Table I — Parameters Used in Fitting /i-Parameters 

Rp, sheet resistance in passive base, 170 ohms/[~l 
Ra, sheet resistance in active base, 6000 ohms/Q 
G E , conductance per unit area of emitter-base junction, 0.00051 mho/mil 2 = 

79 mho/cm 2 
C E , capacitance per unit area of emitter-base junction, 21 pF/mil 2 = 330 /*F/cm 2 
G\, low-frequency limit of ho,, 0.0003 mho 
Cc, capacitance per unit area of collector-base junction, 0.226 pF/mil 2 = 0.0350 

M F/cm 2 
G m , "transconductance" per unit area, 0.058 mho/mil 2 = 8900 mho/cm 2 
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Fig. 5 — Input impedance, h tl or h ic , in ohms. The solid line and circles are experi- 
mental, and the open circles theoretical. Frequency is a parameter along the curve, 
and selected frequencies are given in MHz. 

quency is a parameter along each curve; selected frequencies are indicated 
in MHz. At these frequencies, the calculated /^-parameters are plotted 
as open circles. 

At the higher frequencies, the discrepancy between experimental and 
calculated results becomes significant. There are several reasons why 
this is to be expected. In the first place, the experimental data may be 
in error, since /(-parameter measurements are considerably more difficult 
to obtain at high than at low frequencies. However, it is not necessary 
to invoke this explanation (the traditional one for theorists) for the 
discrepancies to be understood. 

There are two other sources of discrepancy. The first is that the simple 
model used for the transistor is not adequate at high frequencies. 
"Second-order" effects such as transit-time delay in the base and 
collector were ignored; other effects which are small at low frequencies 
could also be important. At high frequencies, the exact shape and size 
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Fig. 6 — Forward current gain, h 2l or h/ e , dimensionless. 
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Fig. 7 — Reverse voltage transfer factor, h 12 or h Te , dimensionless. 



of the transistor is important and the measured transistor may certainly 
differ somewhat in geometry from the nominal. The effect of geometry 
may be seen in the plot for h 22 . The loop in h 22 at high frequencies is 
due to the geometry of the transistor. Calculations done with other 
geometries show that, as the base and emitter stripes are made longer, 
the horizontal width of the loop shrinks; the loop is absent for a transis- 
tor with infinitely long stripes. 

Secondly, at high frequencies the parasitic effects of the header and 
the encapsulation become important and need to be accurately modeled 
if the experimental data are to be matched. Sample calculations with 




0.002 0.004 0.006 0.008 0.010 

Fig. 8 — Output conductance, h 22 or hoc, in mhos. 
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representative values showed that header parasitics had at least as large 
an effect on the /i-parameters as the high-frequency discrepancy. 
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